Cooling trapped atoms in optical resonators 
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We derive an equation for the cooling dynamics of the quantum motion of an atom trapped by 
an external potential inside an optical resonator. This equation has broad validity and allows us to 
identify novel regimes where the motion can be efficiently cooled to the potential ground state. Our 
result shows that the motion is critically affected by quantum correlations induced by the mechanical 
coupling with the resonator, which may lead to selective suppression of certain transitions for the 
appropriate parameters regimes, thereby increasing the cooling efficiency. 

PACS numbers: 32.80.Pj, 32.80.Lg, 42.50.Pq 



Cavity cooling is a recent expression, which stresses 
the role of the mechanical effects of a resonator on the 
atomic and molecular center-of-mass dynamics. Indeed, 
the coupling between resonator and atom gives rise to 
complex dynamics, one aspect of which is the substan- 
tial modification of the atom spectroscopic properties [l[ . 
This property allows one to change the atom scattering 
cross section, thereby affecting, and eventually tailor, the 
mechanical dynamics of the atomic center of mass 
Moreover, the motion of the atom changes the medium 
density, thereby affecting the resonator field itself. Sev- 
eral recent experiments have reported relevant features of 
these complex dynamics 0, 0, HI 0, S S|» Experimental 
demonstrations of atom cooling in resonators [3, EL IE 
have shown, among others, that cavities are a promising 
tool for preparing and controlling cold atomic samples 
of scalable dimensions, which may find relevant applica- 
tions, for instance in quantum information processing Q. 

In this letter we present a study of the quantum dy- 
namics of the center of mass motion of an atomic dipolc, 
which couples to a resonator and to a laser field driving 
it from the side. The system is sketched in Fig. [T] Dif- 
ferently from recent theoretical works on cavity cooling 
of atomic clouds [lOj, here the center of mass is con- 
fined by a tight trap, in a configuration that may cor- 
respond to the experimental situations reported for in- 
stance in 0, EH, El]- Starting from a master equation 
for the quantum variables of dipole, cavity and center-of- 
mass, we derive a closed equation for the center-of-mass 
dynamics. This equation generalizes previous theoretical 
studies [l3|, 14 1 and allows us to identify novel param- 
eter regimes, where cooling can be efficient. Moreover, 
its form permits us to identify the individual scattering 
contributions, thereby gaining insight into the role of the 
various physical parameters. We show that quantum cor- 
relations between atom and resonator may lead to the 
suppression of scattering transitions, thereby enhancing 
the cooling efficiency. 

The starting point is the master equation for an atom 
of mass M whose dipole transition between the ground 
and excited states \g) and |e) couples (quasi-)resonantly 



FIG. 1: An atom is confined by an external potential inside 
an optical resonator, and couples to a laser and a mode of 
the resonator. Energy is dissipated by spontaneous emission, 
at rate 7, and by cavity decay, at rate «. Inset: The dipole 
levels \g) and |e). A and S c are the detunings of dipole and 
cavity from the laser frequency. 



to a laser and the mode of an optical resonator with wave 
vector kh and k c , respectively (|fez,| ~ |fe c | = k). In the 
reference frame of the laser the state \e) is at frequency 
—A and the cavity mode at — S c . The atomic center- 
of-mass is trapped by a harmonic oscillator of frequency 
v along the x-axis, while the motion in the orthogonal 
plane is frozen out. Later on we discuss how the treat- 
ment can be generalized to three dimensional motion. 
The density matrix p of cavity mode and atom internal 
and external degrees of freedom evolves according to the 
master equation dp/dt = Cp where the Liouvillian C de- 
scribes the coherent and incoherent dynamics. External 
and internal degrees of freedom are coupled through the 
spatial gradient of the interaction with the electromag- 
netic field, which scales with the Lamb-Dicke parameter 
r\ = ^Tik 2 /2Mv [lij ]. We assume tight confinement, such 
that the gradient of the e.m. -field over the atomic wave 
packet is small and the Lamb-Dicke regime holds, 77 <C 1, 
and treat the coupling of the center-of-mass with dipole 
and cavity variables in perturbation theory. At second 
order in r\ we decompose C = Cq + X]L\ + i] 2 C2 + o(r/ 3 ) 
where the subscripts indicate the corresponding order in 
the expansion. The term Cq = Cqe + £01 is the sum 
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of the Liouvillian for the external (E) and internal and 
cavity degrees of freedom (7). Here, CqeP = [-ffmec, p]/i?i 
and i7 mcc = hv(n + 1/2) center-of-mass harmonic oscil- 
lator with n number operator for its phononic excita- 
tions; the term C^ip = [Ho,p]/ih + (C s + K,)p, where 
K and C s describe cavity decay at rate k and dipole 
spontaneous emission at rate 7, respectively and 
H = -TiAa^cr - M c a f a + Ti[a\Vl + go) + H.c], with 
(7 = \g)(e\, a,a^ annihilation and creation operators of 
a cavity photon, g coupling constant between cavity and 
dipole at the center of the trap, and laser Rabi fre- 
quency. The first-order correction C\p = [V,p]/ih con- 
tains the operator V = (Vl + V c )(b + tf), which gives 
the mechanical coupling induced by laser and resonator, 
whereby b, b' are the creation and annihilation operators 
of a vibrational quantum, the operators 



Vl = ihtpLQ(a^ — a) 
V c = h(p c g(aa' + a) a) 



(1) 
(2) 



and the coefficients tp^, tp c scale the corresponding recoil 
and depend on the geometry of the setup [17| . The term 
£2 describes mechanical effects associated with the spon- 
taneous emission of a photon, which determine diffusion, 
and state-dependent energy shifts 18,1^]. The modifica- 
tions induced by these shifts in the cavity field dynamics 
have been discussed in [2C}. In the present treatment, 
where 77 1 , they are high-order corrections to the ex- 
ternal potential and hence do not appreciably affect the 
center-of-mass dynamics. 

At second order in the Lamb-Dicke expansion we derive 
a rate equation for the occupations of the number states 
I to) of the center-of-mass harmonic oscillator T3,21|. The 
rate equation is expressed in terms of the heating rate 
= r] 2 (n + 1)A+ and the cooling rate T 

n — m — 1 — 

rfnA- , which describe transitions that change \n) by one 
excitation. Here, 

A± = Re{2Tn{V(C i + i^Vpo} + a 7 Tr{|e)(e| (0o }} 

with T17 trace over the dipole and cavity degrees of free- 
dom and a geometric factor giving the angular dispersion 
of atomic momentum due to the random recoil of spon- 
taneous emission. Their explicit evaluation requires the 
knowledge of the spectrum of the correlation of the opera- 
tor V and of the steady stat e pn of dipole and cavity mode 



at zero order in 77 [18j, [19j, |22|. When the laser weakly 
perturbs the atom-cavity ground state |</>o) = \g, C ), we 



find [21 1 
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where 



5 C + in/2 

HO) 



(4) 



77^ = in 



[S c T v + i/c/2) 



HT") 



Htv) 



/(0)/(=F«/) 



Z K ^ = -n 



g[(k^v + i 1 /2)(5 c + w/2) + ~g 2 ] 

/(o)/(=b0 



with 



f(x) = (x + 5 c + i/c/2) (a; + A + i 7 /2) 



(5) 
(6) 
(7) 
(8) 

(9) 



Equations §5§ and ([I])-© are the main result of this let- 
ter. They have been derived for one-dimensional motion. 
However, since at second order in 77 the three directions 
of oscillation decouple in an anisotropic trap, they can be 
generalized to three-dimensional motion as they hold for 
any geometry of the setup. Their analytic form allows 
one for insight into these complex dynamics, which we 
now discuss. 

The rates are the incoherent sum of three contribu- 
tions. The first term on the right-hand side (rhs) of 
Eq. ©, 7a|7s| 2 , describes absorption of a laser photon 
and spontaneous emission, whereby the change in the 
center-of-mass state is due to the recoil induced by the 
spontaneously emitted photon. This process is shown in 
Fig. [UJa), where the excitation paths are depicted in the 
dressed states basis |±; n) of cavity and atom,which are 
superposition of the states \g, l c ) and |e,0 c ) [2j at fixed 
number of phonons |n). The scattering rate is scaled by 
the geometric factor a and is found after averaging over 
the solid angle of photon emission into free space. This 
contribution is diffusive, as the motion can be scattered 
into a higher or lower vibrational state with probabili- 
ties depending on the overlap integrals [I9j. Diffusion 
can be suppressed, namely term (fj| can become very 
small, in good resonators, K <§C g, and for 6 C = 0. In this 
limit the corresponding dynamics are characterized by 
destructive interference between the two excitation paths 
\<j>o;n) —>■ |±;n) [H, 0, HE]- This situation appears at 
zero order in the mechanical effects and is reminiscent 
of analogous phenomena encountered in the cooling dy- 
namics of trapped multi-level atoms 26] . 

The second term on the rhs of Eq. ([3]), "/{(PlT^' + 
ifcT^'^ | 2 , describes scattering of a laser photon by spon- 
taneous emission where the motion is changed by me- 
chanical coupling to the laser (T^' ) and to the cavity 
(7^ 7,± ) field. Process 77/ scales with the geometric fac- 
tor ipL, which accounts for the recoil due to absorption 
of a laser photon, and is depicted in Fig. [2]^b). The tran- 
sition amplitude T^ <,± 1 shown in Fig. HJc), describes the 
mechanical coupling between the dressed states |±) due 
to the resonator. It thus scales with the geometric fac- 
tor (p c which accounts for the recoil due to interaction 
with the cavity mode. Since the final state of the two 
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FIG. 2: Scattering processes leading to a change of the vibrational number by one phonon. The states \g,O c ;n), |±;n) are 
the cavity-atom dressed states at phonon number n. Processes (a)-(c) describe scattering of a laser photon by spontaneous 
emission. They prevail in good resonators, for c < 9,7. Here, (a) Ts; (b) ?~2 ; (c) T^'^ ■ Processes (d) and (e) describe 
scattering of a laser photon by cavity decay. They prevail in bad resonators, for 7< j,K. Here, (d) T£' ; (e) T*' . See text. 



scattering processes is the same, they interfere 27| . In 
addition, each term is composed by multiple excitation 
paths, and can vanish in some parameter regimes. A 
representative case is found in good resonators, k <C g, 
when the laser is orthogonal to the motion axis (ipL = 0). 
Then, the transition amplitude, T^' + is suppressed for 
S c = v/2, as one can see in Eq. ([7]). This suppression 
arises from destructive interference between the multi- 
ple paths of mechanical coupling via the dressed states 
±; n) and |±;n+ 1), thereby giving a vanishing heating 
rate. Similarly, the cooling transition can be suppressed 
by choosing 8 C — —v/2. 

The third term in the rhs of Eq. ([3]), kI^T^^ + 
iPcT^'^ | 2 , describes scattering of a laser photon by cavity 
decay, where the motion is changed by mechanical cou- 
pling to the laser (T^) and to the cavity (T*^) field. 
The scattered photon is transmitted through the cavity 
mirrors into the external modes, and therefore these two 
processes do not interfere with the ones above discussed. 
They are depicted in Fig. [DJd) and (e) and add up co- 
herently with one another. An interesting situation has 
been discussed in [13j, where it was shown that inter- 
ference between these two terms can lead to suppression 
of the heating transition in bad resonators. The result 
of [IH is recovered from Eqs. ([H])-® in the corresponding 
parameter regime [21| . 

The general dynamics are a competition of all these 
processes. By inspection of the form of Eqs. (J3EJ) , one 
can identify a strategy for obtaining efficient cooling by 
searching for the parameters such that Re{/(z/)} = 0, 
thereby minimizing the denominator of A- , and therefore 
enhancing the cooling rate over the heating rate. This 
condition leads to the equation 



Aopt(^c) 



g 2 + 7^/4 
S c + v 



(10) 



which relates the cavity detuning S c to the atom detuning 
A for fixed couplings and decay rates, such that, when 



A = A opt ((5 c ), the cooling rate A- is maximum. In gen- 
eral, Eq. (fTU)) sets the best parameters for ground state 
cooling, provided that cither 7 <C v or K <C v. Below we 
discuss cooling in the good cavity limit. The bad cavity 
limit will be discussed elsewhere. Figures[3ja) and (b) are 
contour plots, displaying the average number of phonon 
at steady state, (n) — A + /(A- — A + ), and the cooling 
rate, W = rj 2 (A- — A + ), as a function of A and 6 C for a 
good cavity and k <C v <C 7, namely when processes in 
Fig. [DJa)-(c) are dominant. These curves have been ob- 
tained from Eq. ([3]), and their validity has been checked 
by comparison with full numerical simulations [2lj . The 
dashed line indicates the curve A = A opt (<5 c ). It is ev- 
ident that this curve corresponds to the parameter re- 
gion where the lowest temperature is achieved. The final 
limit is set by the ratio k/v, since in this regime k is the 
narrowest linewidth one can achieve for the scattering 
processes [28j |. 

Good parameter regimes are found for relatively small 
values of A and S c , in the range < 5 C < v, A > 0, 
where one has low temperatures and relatively large cool- 
ing rates for a wide range of values. These regions cor- 
respond to novel cooling dynamics, which are character- 
ized by suppression of heating and diffusive transitions 
by quantum correlations induced by the resonator. Here, 
for instance, the potential ground-state occupation can 
reach 99% in a time r ~ 2 msec for an atom in a trap 
with frequency v — 2ir x 500 KHz and r\ = 0.1, whose 
dipole has line width 7 = 2ir x 5 MHz and couple to a 
resonator with g = 2n x 5 MHz and k = 2n x 100 KHz. 
Larger values of k lead to smaller occupations, which 
scale like (n)oo ~ k 2 /v 2 in good resonators for n < v. 
Figure [3] shows also that the dynamics are relatively in- 
sensitive to fluctuations of the parameters in the range 
< 6 C < v, A > 0. This implies that several oscillation 
modes, like the ones of an ion chain [5§|, could be simul- 
taneously cooled. We remark that good cooling regions 
are found in Fig. [3] also for large |A| and S c ~ — v. This 
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FIG. 3: Contour plot of (a) the average phonon number at 
steady state {n)^ and (b) the corresponding cooling rate W 
(units of 10 -3 ^) as a function of S c and A (units of v) in 
the good cavity limit. The darkest region corresponds to the 
smallest values, the lightest region to the largest values. The 
numbers indicate the values at the corresponding contour line. 
The heating regions are not coded and explicitly indicated. 
The dashed line represents Eq. (110[) . Here, r\ — 0.1, = 
ip c — -i/l/2, Q. — v, g = 7u, 7 = Wu, k — O.Olu. 



is the so-called cavity sideband cooling regime, discussed 
in 14|. Here, enhancement of the cooling rate is achieved 



by means of resonance with the very narrow transition, 
at large |A|, between the dressed states of the cavity- 
atom coupling. This regime exhibits a critical sensitivity 
to variations of the value of S c with respect to A. 

To conclude, we have presented the explicit form of 
cooling and heating rates of the center-of-mass motion 
of an atom, trapped in an optical resonator and driven 
by a laser. This result permits us to identify the basic 
processes determining the center-of-mass dynamics and 
the parameter regions where cooling can be efficiently 
obtained. Novel regimes have been identified, where the 
motion is cooled by exploiting destructive interference in 
the multiple paths of the mechanical excitations. The 
resulting interference bases itself on the discreteness of 
the vibrational spectrum, which is the same for the dipo- 
lar ground and the excited state. Some analogies with 
interference _processes in cooling of trapped multilevel 



atoms 
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are due to a similar dressed state struc- 
ture at low excitation, nevertheless the mechanical ef- 
fects induced by the resonator give rise to a wider wealth 
of phenomena. Dynamics will be substantially modified 
when the external potential depends appreciably on the 
internal state [20L I22I ]. In future works we will investi- 
gate how the dynamics are modified by collective effects, 
when many atoms are trapped inside the cavity. 

Beyond applications to cooling, these results allow one 
for gaining insight in the complex dynamics of the me- 
chanical effects of optical resonators on atoms, and their 
wealth of phenomena could be exploited for implement- 



ing coherent control of this kind of systems. 
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